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I. INTRODUCTION 


The Short Range Air Defense Command and Control (SHORAD C2) system con- 
sists of the equipment necessary to provide C2 capability for Army divisional 
and non-divisional Air Defense Artillery (ADA). It will provide C2 for all 
echelons of a SHORAD battalion. The network will consist of displays, data 
processors, interface equipment, and sensors. The ssjor objectives of the 
system are: : 


- Provide automated alert to all SHORAD battalion elements. 

- Provide for automatic exchange of command information between 
SHORAD battalion elements. 

Provide weapon cueing to selected SHORAD battalion elements. 
Provide track information from the various sensor inputs to 
SHORAD C2. This involves track management functions which 
include track file management and track smoothing operations. 


Track files from external sensors such as PATRIOT, TSQ-73, HAWK, ADEWS, 
and AWACS feed data into high level modes of SHORAD c3, Track management func- 
tions include the maintenance, correlation, triangulation and smoothing of track 
data. This report descibes a technique for smoothing noisy track data using 
a GHK discrete time-recursive state estimator program. The GHK estimator is 
useful in smoothing noisy track data obtained from noisy measurement sensors 
in applications such as the c2 system for SHORAD. This computer program was 
developed as part of an analysis of SHORAD track management functions. 


A detailed description of inputs and calculations will be included in the 
discussion to assist the user in choosing the original values. In addition, 
examples of inputs, resulting calculations, and plots will be provided in the 
results to give the user a thorough understanding of the output. Also, the 
mecessary system assignment statements to aid the user in obtaining a hard 
copy of the results will be included in Appendixes A - C. 


TI. DISCUSSION 


Initial values are input to determine the initial position, velocity, and 
acceleration of the target. The program then estimates the current position, 
velocity, and acceleration at each discrete time, based on initial conditions 
and a predefined predictor/estimator logic. Noise and various maneuvers may 
be added to the true position to determine how accurately the estimator 
handles noisy sensor data on maneuvering targets. Plots of true and estimated 
positions may be obtained to analyze the adequacy of filtering logic. 


NOTE: Name of output entered is used as graph title. 
A. Inputs 
1. Main Program Inputs 
VARIABLE FORMAT DESCRIPTION 
SEED 15 Used to generate random numers. Must be an odd 


integer 0 < SEED < 32768. The SEED is used to 
generate a unique set of random numbers. 


™ 
TMAX 


¥6.2 
E12.5 


' E12.5 


E16.9 


£12.5 


¥6.2 


F6.2 
E12.5 


£12.58 


F6.2 
F6.2 
F6.2 
13 


Estimate of position. 

Time derivative estimate. 

Eetimate of acceleration. 

BETA + BETA may be entered from the keyboard or 
computed. Please refer to calculation section 
for computation. 

NOTE ~ G, H, and K are a function of BETA. 
Please refer to calculation section for 
computations. 

Noise SIGMA. 

Noise MEAN. 

Initial true position. 

Initial time derivative. 

Initial time. 

Time increment. 


Time of maneuver start. 


Number of computation steps. 


All inputs except IMAX and SEED are floating point. 


NOTE: An option is available so that if one chooses to rerun the program, 
he will not have to reenter the original values. When this option is chosen, 
the user can change only the specific value(s) of interest to him. 


2. Inputs for Subroutine BETA 


Te 
DELT 


B. Computations 


F5.2 
F5.2 


Time constant 
Time increment 


1. Subroutine BETA~BETA is computed as follows: 


E = 2.718281828 


LAM = 1./TC 


X = LAM x DELT 


BETA (B) = E**X 


so BETA = .967216134. 


2. 


Subroutine GHK 


Gs 1. = Ba83 : 
H = 3./2. x ((l.- BE#2) x (1. = Y)); (Y = noisy position measurement) 
K 2 1./2. x ((1. - Y) $#3)) 


3. 


Main Program 
a. Target Computation 
Target acceleration is computed as follows: 


ZDD = 78.5 (T-TM) 


User inputs 
T - Initial start time 
T™ - Time when maneuver starts 
Target velocity is modeled by: 
Z1 - Initial true position 
ZD = Time derivative 
T - Start time 
Ir (Tt GE. TM) Zz 2 Z + 0.5 x ZDD x (T=TM) #82 
(0.5 x ZDD x (T-TM)##2) = Position change due to maneuver 
be Noise Modeling 
Subroutine GAUSS - Generates random numbers 
GAUSS (NO, SIG, RMEAN, WG) 
NO = Seed 
SIG - Standard deviation 
RMEAN = Desired mean of the normal distribution 
WG = Value of the computed normal random variable 


Y s Z + WO (where Y = noisy position veasurement) 


C. Logic 
1. Estimating Logic 


XBAR = XHAT + G x (Y-XHAT) 
XDBAR = XDHAT + (H/DELT) x (Y-XHAT) 
XDDBAR = ZDDHAT + ((2 x K)/(DELT**2)) x (Y=-XHAT) 


2. Predicting Logic 
R = (DELT**2)/2. 


XHAT = XBAR + XDBAR x DELT + XDDBAR x R 
XDHAT = XDBAR + XDDBAR x DELT 
XDDHAT = XDDBAR 


IIL. RESULTS 


Each value entered is repeated to provide the user with an opportunity 
to change or correct value(s). Moreover, one has the choice of printing 
‘detailed results for each step (Table 1), or a tabular summation of 
results (Table 2.) 


Examples of inputs and corresponding results are shown in the figures 
that follow. For odd numbered figures, the time constant is 3. and DELT is 
el, so BETA = .967216134. 


IV. SUMMARY AND CONCLUSION 


This system shows how accurately a target can be tracked given the 
appropriate input values. Therefore, it is very important to use the initial 
values that will yield feasible results. Any change in the main program input 
values on pages 1 and 2 will change the results. 


When BETA = .8 (even numbered figures), the program is able to track the 
displacement vectors fairly accurately. Yor example, in Pigure 7, (BETA = 
-967) there was more lag in the result, but the curve was smoother. In Figure 
8 (BETA = .8) it tracked more accurately, but there was more noise. The 
program was able to track the velocity vectors much better with a BETA of .967 
as opposed to a BETA of .8. The results depend on how closely the initial 
filter states correspond to the true state values and on the filter damping 
design, which is determined by the choice of BETA. 
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MAIN PROGRAM 


THIS [S A Gen EST +44T QW 

OIMENST Or ZaCy tty, weaPe (ED, Zoe hb teed, aPtAF ACH.) 
OLMENSTUN ZOMeRC Pee TT Chet), adraddde.) 

CHAZACTES way NAME 

REAL “,*V 

INTEGER C1 

R zs (DELTe#2)/2, 2 

w Zu. 

ci ze 

oO zs 1 


FORMAT( 1s Sane INOUT SECT ION°) 


FORMAT Le o eee eee eR Aka RER HERRERO ERAAERHTHHEERIRERHEER 


Leerareranheanenteexnareeenan ReRReHEeReRteeEPeKEKRatAuvananae’) 


tS 


TYPE 30 

FORMAT(° 1°, Seeemee ALL YALUES EN FERED ARE FLOATING ERENT ee ks 
OTHERWISE LM ICaATENxawre’s/) 

TvPE av 

FoRATC LS, “PLEASE ENTER GAME TH CaLL YRIS SET *8 Curt *s) 
PEANCS, Sut) sae 

FOR AT (AA 1) 

HALT. T ot, sae 

BOR aT( LAK, *vaANME JF LVUTPUT LS %, 48) 


TYPE 8H 

FOR“ATC( UX, °SEED “UST SE Ati 109 INTEGER %//, 

© PLEASE ENTER SEED suttEr FOE SUBROUTINES) 

REAC(S, 990ER® = Fu) cid 

FORMAT (15) 

PRIUT La, 40 

FORMAT (2%,°SEED NUMMER [5 °, 157) 

For4atT (9X,° SEED FOS TerI§S RUN [§ %, 15/7) 

To SEE IF SEED IS wITHIN SPECIFIED BOUNDS 
TFOND «LE. 4 .GP. NOU 4uT. 327608) THEN 
PRINT €, °ERROQReweees SEED HUST SE HETAEES + 15> 327HA" 
Go fT) 7a 

ELSE 
COuTINUE 

ENDIF 


C CrECK TO SEE IF SEED 1S "POD 


tev 


wet 2 ase 

wey Bw Nos2!, 

TFECt ry Ewe mt) Tees 
PRINT 8, eee FORI) neon” 
SU TT) 7A 

ELSE 
CONTINGE 

ESOLF 

WATTECOe1-) 

wWRITECOs 20) 

WRITE (oe64) NAHE 

aPITECO,119) Oa 


IF ("8,E9.1) GO Te 175 
TYPE 12.) 
FORPMATC( 1X4 CQAPBESTINATE SF SOSTTION° SS 1X0" aNvee 2 TE F CEs rvaly] 


LVE ESTIAATE'/, 14%, % ADOQAWSESTIMATE OF ACCELE@ ATI ur °7) 


63 


27y 


TYPE tau 

FOR AATC IX, °PLEASE EXTER £oae, CROAK, BUD RUN sAR%) 
Cis. : 
REAr(Se 1S ',ERRFSeL Se) sRAR, YOR AY, XNDRAY 
FORMAT(FAL2,2(E12.5)) 

FQRMATCIiX,? KGaRs %,FAL2,% wUnbes %,E12.5,° Kedears %,E16.5/) 
FQRMATC LOS, % KRARS %fFAL2,° KOOKS 9S E1205," WyRses %, 
€12.5/) 

TYPE 163, XBAP, XURAP, XKDDBAR 

TYPE 194 

PORMATCIXG “ENTER 1 TO CraiGe VALUES JUST ENTE ED?) 

accePr 24%, C1 

FoRaT(It) 

IF(CL EQ. 1) GO T9 134 

WRITE COs 72) HAI, XOBAR, KUDEAR 


MCLE VALVES 
XBV 3 X4AR 
XDEV = XDRAF 
aCORV = XUDGSAR 


IF (NR .Ets. 1) THE 
FQINT 212, Ge tty A 
PRINT 215, + 
FORVYAT(2xk,°S8ETA TS %sEtee 2/7) 
Eye AT (LOA, GS Moe L dee”? AS “FEL tele” KE CoC tdL7/) 
PRINT @&, “ENTER 1 TC CHa.GE FETA? 
aCCEPT au, Ct 
IF(C1 eE9. 1) CALL Grn(Ge oe Ke 5) 
CONTINUE 
ELSE 
CALL GH (Gy Me Ke 3) 
WRITE Cos2ti*) Ge Hy K 
ENOIF 
MOLD VALUC CF BETA Gy Hy, AND K 
Gv 3s G 
hy 
KV 
AV 


ro) 
« 
3 


IF (82 .E4%, 1) GO Tu 255 

TYPE 23 

Cis ” 

FORMATCIX. ENTER GAUSTAN “NOISE SIGMA 4%) NOISE “F4AN%) 
PESO (5,204, EF 2522::) SIG, IMEAN 

FORMATCIXe “NOLSE FACTIM IS %,Fo.2e” YEN! IS “ep Fn.c/) 
FOR ATC Lor, °MOTSE FaCTinh [5 "F602," “EAN TS %,ERL Os) 
PAlT:T 24d, SIG, RYEAS 

FORMAT (2FO.2) 

TYPE 193 

acCEPT ats, 0) 

IFCC1 o£556 1) GO TO 225 

aA{TE (0-25) SIG, F4Ear: 


SIGv 2 Sf 
PMNY & RME AR 


IFC UR Ewe 1) 69 TO 325 
TYPE 27V' 
FOR MATOCIXs “ZISINITIAG THUE POSITION? /IX, °ZRLSTI ME LEST VAT IVE ls / 


1 xe °ZUOSACCELERATION?7) 


28? ci = 7 
Pal u TYPE 38) 
30a FORMATC IX, PLEASE EXTER Zi, ZO21, sNu 2ML"%) 
QEAV(Se See ERPs294) Z1, Zulse 209 
31¢ FOR“ATCIN,? 212°, EL205e"% ZOIE %-E12656% ZF %e 212.57) 
32a FORMATCL5SX,*% Zia’ cEl2e5,"% ZO1S %e61205,% ZS °,F 12.57) 
325 PRINT 314, 21, ZD1, 200 
TYPE 19% 


ACCEPT 200, Ci 
TF(C1 E96 1) GO TO 254 
WRITE € 6¢32@) Z1, Z01, 250 


c 
Z71V¥ 3 Zi 
ZO1v 2 ZD1 
Zo00V = ZNb 
c 
IF(NR EG. 1) GO TO 445 
TYPE 331 
332 FORMATC I, °T = TIME%, /1x,7DELT = TEME [T° CKEMENT?, stare T* = TY FF 
IMANEUVER%//) 
Cc 
Cc 
TYPE Bac 
340 FORMATCIX, “IMAXSNUMHER OF COMPUTATIONAL STEPS°%/) 
380 TYPE 36:! 
ci = 
360 FORMAT(IX, °PLEASE ENTER T, MELT, TY, AML JAY CP. TEGEF)%) 
. READ(S,375,ERR2357) T, DELT, TT’, IMAX 
37e FOR“4AT(3F0.2,13) 
382 FORMATCIXs *IMAXsS %, 13/) 
398 FORMAT C LOX, °IMAX= %¢ I3/) 
4ee FOP“ATCIX,% T= °oFaeds*® VELT=S %,F5.2,° TN Ae 2s) 
ated FORMATCISAS® Te "¢Foede”® DELTS °oF beer” Tie %e Fb. 27) 
ais PRINT 33, [MAX 
POINT dv, T, DELT, TT 


TYPE 19s! 

ACCEPT ane Ct 

TFCC1 E056 3) GN TD 35V 
ARITE (Coe39N) [MAX 

wQITE Coe4i*) Tye, DELT, Th 
AR ITE (60253) 


C HOLD VALUES OF T, CELT, Tt, At [vax 
Ty sf 
Tuy s TT"! 
CELVY = JELT 
TUaAKv = [tax 
c 
HR TTE Coe 427) 
uev Foxear( 7173 
c 
C xaaaeaaaeaeraner COMPUCTAT[UN SECTS 1% aeeeanerenetentens 
430 TYPE 435 
Cycece 
435 FORMAT Ctx,’ ENTER 1 TU PRT IT RESULTS FOR EACH STE %) 
READS, 2A5,ERR 243.1) C1 
c 


‘ CALL GAUSS (NO,SIG, OME At, 015) 


CQO tae Tat, ray 
f= Tf ¢ JELT 
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c : 
Cc STOWE VOLES To 8€ PPINTEI, Tre at ADLAY 
c 
ZHCT) = 7 
MBA CT) = X3AQ 
ZNHCT) = 29 
XYOBAPHC IDS De aR 
ZO04¢i) = 200 
xNOem (1) = xOngar 
TT(L) = fT 
c 
4ayu COMTI Ve 
Cc 
c eeeneanwes PRIGT VALUES STORED Ti ARYAY aeanneeaneee 
RITECO, S30) 
Sau FORMATO TL ESTER am Zh Le ey RAE LK LOL a Oe, 
1 Axe "ZOD g hix,y PXONGARTS/) 
Cc 
oe $9: JP = t,T41Aax 
MPITECO, of) JP, ZAC IP), Kotte IP), ZH IE), UC erN-( SF), 
t ZOO IC JP, KARR IP) 
ota FOR ATC UXe Ty, 1X, E1265) /7/) 
590 Contluue 
c 
C CmECK TU SEE wrETNES OLOTS ALE wAtTe™ 
fpe* 
POT. 0, SLGES LTH PUT vata? 
ACCEPT ec, C1 
IF(C}Y 2&9. 1) THES 
Cc CALL PLOT PROGRAM 
C CHECK TO SEE ehETHER TI PLUT DISPLACE RST HR VEL ICTY vey toes 
603 TYPE 62S 
OAS FuRmat (% ENTER 1 TD PLOP Ve TSCLACEME Dh VEE ty 
1° FUTER 2 TO PLOT VELOCITY VECT wes, 
1° ENTE“ 3 TO PLOT ACCELERAT I. we CT GH Ys) 
Clas 
ACCEPT ay, Cl 
IF (Cy eB 1 oad, C1 les é oA, C4 . Fe 3) Tess 
PLDT #, “PLESSE ENTES 13 7R Aun 3° 
G) TO @.'3 
ELSE 
Co aT Tele 
Eile 
IF(cCt .€&, 1) CALL PEST (TI, reaus, 2a, trae, sch) 
TFCCi JE. 2) CALL PAST CTT, st bane 2 myo er, so 2) 
TFCCL Ff, 8) CALl Fett CTT, Lit aya ryt ocr, eG} 
Cc Ew CE SUT wher pre 
ELSE 
COUN EAs 
€ OIF 
c 
PRI:T mpuy, (ACE 
WETTECO,ALS) Ath 
61{@ FIWMATCL Xe SENG OF Jak Fok Corp ae "ae ') 
c 
Cc CHECm Flite ASAT HER wert, 
TYPE 62" 
‘he oS | 
et Fem «atl C*° EF ute® { Th thw h AT TERS we | 
ACCEPT eu, xb 
TF (rx SE eo. 1) Tr&e 
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a 


ary 


7SINIT 


TF (N32 

TALIZ: 
Xia 
X93 
xDD 
$ic 


eae 


Vays 


gos 


&G = 


3a 


ee ar 


REFERENCES 


“"Bistatic Tracking and Trajectory, Estimating Techniques," 
Henry C. Holman, Dissertation, Southeastern Institute of 
Technology, 1981. 


"VAX 11/780" User's Guide. 
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APPENDIX A 


OPERATING PROCEDURES 
1. Procedure for running program: 


a. LOG ON TO VAX 780 

b. TYPE "ASSIGN EST.OUT FOR 006" 
c. TYPE "ASSIGN EST.TEX FOROOL" 
d. TYPE "RUN EST" 

e. FOLLOW PROMPTS 


2. Procedure for obtaining hardcopies: 


a. TYPE "PR EST.OUT" 

b. OBTAIN HARDCOPY FROM OPERATOR 

c. LOGIN ON TEKTRONIX NUMBER 

d. TYPE "EST TEX" 

e. PROGRAM SUPPORT CONTACT - SANDRA BRAZELTON/876~2295 


A-1/(A-2 blank) 


APPENDIX B 


LOGIC BLOCK DIAGRAM 


B-1/(B-2 blank) 


WVYOVIG ADO1d DLOO'! 


APPENDIX C 


SOURCE LISTING 


C-1/(%-2 blank) 


30 


35 
49 


590 


49 
¢ 
¢ 
70 
80 


99 


106 
110 
c CHuECK 


a 
is) 


wars BQ9G2a™ 
THIS 23 3 GBR ESTIAATSR 

DIMENSION ZHC16C), XehaHn(199), LONCTID), KIZAPRC1I9) 
QIMENSION ZOONCTIOI- TCI), XISEHC1OO) 

CHAPACTER #4C NAME 

REAL Ks KV 

INTESER C1 

as Q. 

cra C 

nO = 1 


FORMATC* 1°, 3GXe° INPUT SECTIONS) 


FORMAT (1X s * ea arerarerrrerretaeeereeatenketharennNakTerihy 


1 ee ee eee ReERHOARHERAAHTARTEKERHREHeEHROeReECeTReETERDKUAREteR!) 


1S 


Types 30 
FORMATC' 1°, *eawewe ALL VALUES ENTERED A2E FLOATING POINT --UNLES 
OTHERWISE INDICATE Deere" //) 

TyeWs 4G 

FORMATCIXs*OLEASE ENTER NA¥E TI CALL THES SET OF CUTPUT*,) 
READS 4 59) NAME 
FORMATCASQ) 

PRINT S50, NAWE 
FIRMATCTIXs NAME CS 


mT) 
a 
Cc 
“+ 
uv 
c 
4 
_ 
“ 
- 
. 
2 
a 
=) 
~ 


TYPE 3c 
FORMATCOXs°*SEED AUST FF 
" PLEASE ENTZA SEZ NUME 
READCSEZ,IOFZEW ZB 70) NO 
FORMAT (15) 
PRINT 1Jde NO 
FORMAT C2Xe*SEED NUMSER 75 *, 28s) 
FORMAT (39X,* SEED FOR TRIS PUN TS *FI5//) 
TO SE IF SEED IS aITHI®S SPECIFIED BOUNDS 
IFCHO eLEe JG Ste lS 0 Te FC75R) THEN 
PRINT we SRROR|mMMM N= SEED MUST 32S BETWEEN C AND 32768° 
$0 TS 735 
SLSE 
CONTINUE 
SNOTF 


AN 000 INTEGER s/, 
6? FIP SUSROUTINE*/) 


Chec< To SEZ TF SEED IS Jod 


tv 


S ESTIMATE so 1Xe* 43595482 


NH1 = NC/2 

ANT = NO/2. 

TECNHT 2it%e YNT) THEN 
PRINT tp PemnZRAGIMKMt 


20 To 73 


CONTINUE 
ENOIF 
@R®ITEC6219) 
@RITE (6629) 
WILTE (4540) NAME 
eRITECe119) NC 


ITF (NF.Z%.1) GC TO 175 
TYPE 12° 
FORMATCIXS* KSA zz 


TivaTZ OF SCTITIIN'/,1K," XDFAR = TIME JEVIVE” 
ZESTIMATE JF ACCELERATIONS) 


o 


215 
219 


2793 


TYPE 1469 

FORMATCTAs PLEASE ZNTER AEAT/NDIAI, AND XODFAR*) 
c1 = C 

RENCE S15 I-/ERRZ139) KBAP, XOBAI, KODBAR 
FORMATCFZ ocr 2CE1E05)) 

FORMATCVAs* XBARS 'oFSe2e" XIBARH "FE12.Se* KOOEAI™ °F 2172.54) 
FORMATC14N7" XSARE "oFSecs" XOGAPS %5512.5e" XOCAARS %, 
£12.5/) 

TYPE 14504 X9A7, AXDSAR, XIOBAR 

TYPE 130 


FORMATCIXS* ENTER 1 79 CHANGE VALUES JUST ENTERED*) 
ACCEPT 209% C1 

FORMATCI1) 

IFCC1 223-2 1) GO TS 1735 

WRITE (47179) XEAR, KDZA5- AGO SRR 


WOLD VALUES 
X3V0 2 X3AR 
XDcV 2 XOBAR 
Xoosv = x0094" 


IF (N® 229. 1) THE 
PSINT 2104 Sse Her © 
PSINT 218% & 


FOTMATCIN,*3ETA IS 87215.°/) 
FORMAT (ASX, "Gt S,Elee7e" H® Mp itee7o® KE *FF74.7/) 
D9INT #, “ENTER 1 TS CHANGE SETA 
accePT 206, C1 
YECC1 22Q. 1) CALL GH<(Gs Hse Ke 3) 
CONTINUE 
ELSE 
CALL GHK (Gs We Ks &) 
@RITE (6rc10) Gr tho X 
INOIF 
HOLD 
SV z 
MV 


ALUE OF SETA SG, He AND X 


v 
3 
] 
< 
S$ 


Yew 
dea 
ave 


FIRMATCIX,ZTENTES SAUSTAN N 
2202 (5,250,299 S625) SiGe 
FOPMATCTASTNCISE FACTOR T3 
FIAMAT CTX, NCIES Factor T5 
POINT 2ede sieve F4isii 
FIRMATCTFE 5.2) 

TYPE 199 

accer?r 20C, C1 

TECCT 222. 1) GO TI 229 
@aITE (62250) SIG, FIMZAN 


c 


SIGMA AND NOISE YEA) 


ua 


SiGV 2 32% 
aMNV = PMEAN 


TFONR eSae 12 SS TS 325 
TYP 276 
FURMATCIXS LIPINITIAL TINUE POSITION TX, *2O1aTIMe DERIVATIVE’ os 


Cm 


360 


373 
363 
390 
209 
410 
415 


TVetIOOZACTELEAATICON' SD 

c1 *« ° 

TYPE 700 

FIRYMATCIA,*PLEaSs SITES 241, 201, 32ND 29D*) 
PEANC$4150,29A226G) Zt, IDI, 298 . 
FORVATOC1AS® D1z*,E12oSe* LI1s %,212.55% 2908 %4F12.5/) 
FORMATC15Xe? 21e8 set 2e5e* OTB "212087" 2998 *E12.50) 
PRINT 31G, Z146 314 209 

TYPS 136 

ACCEPT 2004 €1 . 

IFC(C1 2&8. 1) GO TO 2730 

wAIT2 € 64320) 21, 291, 7200 


itv 2 21 
zD1V 2 251 
weov = 200 


SFQNR 283. 1) GO TI 415 
TYPE 350 
ECRYATCVAS OT = TIMEN,JIXSTOELT = TIVE INCSEMENT' A /1Ket TM = OTE 


< 
wm 


TMANEUVEITS/) 


TYPE 769 
FORMATCVXK, IMAXSNU4E TI OF 
TYP= 350 

ct 2 ¢ 

FOR“ATCIXs PLEASE ENTER Te 
REAOCS »370QeERAz750) Tr SEL 
FORMATC3F 6. 2,12) 
FORWATCIXs *I¥AKE 8550) 
FOPMATCIGXs *IMAKS $4237) 
COAMATOCILS' Te "eh aede* JELT= "oF 502" TM %,F402/) 
FORMATLYSXe* Te “st aete* SELTS 8e85e2e9 THR NF FSe2/) 
Palyy 330, Imax 

POINT 4594 Te JELT, T¥ 

TYPE 199 

acCEPT 229,C1 


QMOTATISNAL STz9S5"/) 


ao 


SELT, TM, ANC T4AX CINTESES)®) 
o TM, TMAX 


-? 


. TECCT 263. 12 GO TI 358) 


APITE C5377) IMAX 
aQAITE (4,410) Te Jelts T¥ 
a@lTe(se.c)) 


C AOLD VOLJES TF Tr DELT+ Tse AND TAY 


Ty = f 

Tv =z 1% 
SEuv = ELT 
2MAKV = T¥4XK 


C 
eFlTilses2es) 

6e5 sevarc*t*) 

¢ 

C wenenenenernaree COMPUTATICN SECTION wrewweweecereneren 

430 TYPE 435 
c1 2 C 

@35 FIAMAT C1XF" ENTE? 1 TI PSINT RESULTS FO? ZACH STEP °) 
REANCS S ZG2/ 72868455) CI 

C 


O99 640 Tal, I>ax 
tT 2 # SELT 
IE CF ek Te T¥) 22 TD 657 


co5 = 72.8 © (TT) 
TF C7OS 2 5E. 4561) 259 = 4,1 


250 7) = 251 * 200 & (TT) 
zZ = 24 © ZOeT 
LeCT oGte T4) Z 2 Le Set * LOD & (THT) ee? 
¢ 
CALL GAUSS (NC OSTGsP MEANS dG) 
Y 2 2 + 4 
¢ 
¢ 


KF CI 2E£3- 1:39 GO TE 450 
Ceweeenawewe ESTIMATE AT CUIAENT TIME eevereeneszer 
COREA AETHRAARPEHRARETRETRTERTEO NATTA HECTOR LERTE KER OeP ART D 
C NEXT THREE STATEMENTS ARE SKIPPED CN FIRET PASS od 
C ee RRM OHARA TAT AETEHROTRERTE RA ERERERRERRD ETAT eH eee 
XSAR = XNAT * 3 & CY¥“XHAT) 
MOZAR = KOHAT © CH/TELT) ©& CY“XHAT) 
KODSAR = XOOHAT # CC 2 © C)SCDELT##2)) © CY—XHAT) 


Crewwerevewevenenrene ROReTeTewzwenaeerterereerevareeneverreretankee 


¢ THIS £8 THE LOSTC FOR TRE ELTAST PASS « 
Cer ee Reet eRe TATA EN EAHAHRERERA RE THTHEH TETRA PENRAED OWT ERANT EH Ee 
C 


Q 2 (CBLTewt) s/c. 

66d KYAT = KIAF © ADSAR © O 
KIdAT = KIBAR © ADQIAP «& 8 
XODHAT= XODFAR 


LT * xKoDoHe « 9 
LT 


¢ 

Creaeaneeereeecaeneee PUTHNT QOUTINE stwrewtetvenrarewenrne 

C 

C ewew CHECK T7 SEE IF ALL RESULTS AME TO 2E POINTED wen 
¢ 


IF (C1 224. 0) GO TO 47C 


w 

CC wewe SKIP THIS SECTION CF C1 = GO waeeewee 

\ 
WRITE (46 B3C)2 

639 LFORMATCTIN,VESULTS FOR STEP ',i3s) 
ePlTTE C4, 69C) 

£96 FORMAT CTX eT 8 V3ko* KIAD 86 5X4" XDOSAQ §,5Ke" KODAAR") 
@P2LTS C9, FIDIVT, KZA, YDEAI, KOTIAS 

§ac FIFVATCV ADE Teco Dre 1 2 eso QhoeVeeS ether il2054/) 
e7ITiCs, $19) 

513 FSRMATCSKe* KMAT "59K," XOHAT 6, 35X4" KOSHATS) 
e22TECa, SLT) KYAT, KIMATS KIOMET 

§2¢0 FIRMATOCS ASS CEV2 0305000) 
e2sTitss, $525) 

53) FIIMATCEAS 2 MLAS, lot, pb ke 290") 
ePlTidss, $25) %t4 236 2d0 

$e FIOMAT COX EV ee So eXrilinirekoil2e5ss) 
wesrTé(es, $523) 

559 FORMAT C3XS*XBAR = 2 8, 9K ep LDSAANTO OTK eo XODRAINZ00®) 
APSTECS55459) (Xr aAMZ) A CKDSAINHIOIZCXOOBAR ~ 190) 

£60 FIVIATCEKA 3CE12050 2X00) 


aPITI(er577) Yo CY¥oUnAT) 
73 FIAMATCT as Yt Ss sl ieis ite CYR KAAT) © *FET2.5//) 


5 
c 
C weweerne END DF SECTION TO 32 SKIDOED weenee 
c 
c 


aweearwene 


CHECK TO 


— 


CONTINUSICIT 2 © 


STO®E VALUES TO :E& 


ZHCT) 2 
XZAQNC7) & 
ZOHCZ) = 
XDEARHC(CI)= 
ZTOOHCI)D #8 
XOOSHC(I) # 
TTCI) 2 T 


CONTINUE 


aAlTECs, 


PRIMT WALUES STORED IM ageay 
$3 
FORMAT CPN1%, 


30 CNLY FPILNT VALUES STCRED IN SAPAY 


PIINTES Ty ane ARPAY 
z 

X2Ak 

zo 

XDSAR 

790 

XOOBAR - 


aReeReCRR Ree AES 
3) 
STEP har Zt y Oke  XIAR oTOKs ZOE re TOXe  XDSAR', 


eXe*ZIO 10K *XOOIAR SSS) 


do 539C J? Vel *aXx 
eALTEC4Sse 452) 58, LHC SPD KPAPHCSPIDV, ZOHCIPI, NTSEWMCIP), 
LOOHCIP), KDOSK (4?) 
FOIMATCING E3030 Ke it 2.5/1) 
CONTINUE 
Sz eWHETREF PLITS AGE WANTED 
C129 
POINT w,* ENTES 4 TS SLOT JATA® 


ACCEPT 200. C1 


FCI 


eFde 


1) THEN 


CALL PLOT F90GIA4 


CHECK TO SEE AHETRER TO PLST OTSPLACEMENT CR VELOCITY VECTORS 


TYPé 405 
FORMAT CC® ENTER 1 TO PLOT OLISPLACEWENT VECTOR®/ 
ENTER 2 TO PLOT VELOCITY VECTOR’, 
ENTER 2 TO PLOT ACCELEVATION VECTORS) 
c1 2 9) 
ACCEPT 2C2, C1 
LECC1 New 1 eANDe CT oNEe 2 eANDe C1 oNEw 3} THEN 
PPINT #, "PLZaze ENTER 1:59 2 OP Tt! 
$c To $63 
EL3? 
LONTINGE 
ZNOIF 
TFCC1 22%. 1) CALL PEST CTT/CESEM, THe TMAX/NAME) 
TECST «E52 7) CALL CEST CTT ZAOSARM, TOM, I MAK, NAME) 
IECC1T eise 5) CALL PEST CTT eT DOW se K00 FH TV AK Se NAME) 
af 21 Av FQUTINE 
hae 
CONTINUE 
SUOTE 


POINT o1ce 


NAME 


WOLTE(S,6195) NAME 


FORMATCIXS END OF J33 


FOR PROGRAM °,430) 


CHECK FS2 ANCTHER AUN 


TYyPz 429 
Na 3 3 
FORMAT 
ACCEPT 


(a 


entes 1 TS 
2950, NR 


“AKS ANCTHER PUN °) 


i 


sF ¢ 
JETNITIALT 


N9 2.236 1) Tre 
fe NAMES TS Sf 
XSAR 2 XaVv 
MOVAR = XDEV 
ROOSAQI = XEOSV 
3tG 2 SIGV 
RVSAN = VENV 
21 2 Z21V 

Zot = ZDtVv 

z85 = LooV 

2 av 


~ANA xX OOD 
wu unt 
x 
< 


TY = TMV 
OELT @ OELYV 
2MAX 2 IMAXV 
Ge Ta 5 


CONTINUE 


c-8 


~ ee eee we ee 


ace — ee ne a a 


94 


SUBROUTINE PEST 


SUBROUTINE PEST CTT+ X01, XN2, IMAX, NAVE) 


C THIS SUSROUTINE PLOTS TIME VERSES TRUE POSITION 


CHARACTER #828 NAME 
CALL INITT(120) 

CALL BINITT 

CALL NPTS CIMAX) 
CALL CHECK (TT,X01) 
CALL OSPLAY (TT,Xx91) 
CALL CPLOT (TT,XD2) 
CALL LABLIT (389,NAME,31%,704,1) 
Call FIVITT (4, 71.9) 
RETUR’! 

ENO 


SUBROUTINE BETA 


Sj82%0uTIme @eta tH) 

THIS PP IGHAM CALCULATES Jets 

REAL LA” 

PRIT x, (€ % SLEASE c' TER TIVE CA.STANT %) 
wmeA™ ( 5,0.,E"% = 1°) TSE 

LASS = 1/7 TC 

PRI-T 4%, Lan 


Ee aATCix, ? THE Valet IF Lasara IS %,Fo.2) 
PwP ot we, (C% PLEASE F UTES VFL_T %) 

ta) (Spo*,ka* = F:) urUT 

Fo? At € FS, 2) 

Pe a ee Ge et 

Bom ATC Lee Tan val. 2 UF CELT IS % F522) 
FS 2.71724), 4% 28 

% 2 L444 @ ELT 

BUTT Say. 

FORMATOCL 4, °LAMSUA #DELT = %,&16.9) 

os ©& ee Kk 

RETUP 

El 


NAnNanNnaANaNnaNaNNaANANNAaAaNNnaANANaANnaAanaaNnnNnaNaaNANNANAnNAn 


5@ 


SUBROUTINE GAUSS 


SUBROUTINE GAUSS 


PURPOSE 
COMPUTES A NORMALLY DISTRIBUTED RANDOM NUMAER #ITH & GIVE 
MEAN ANO STANDARD DEVIATION 


USAGE 
CALL GAUSS(IX,SeAM,V) 


OESCRIPTION OF PARAMETERS 

IX ©IX MUST CONTAIN AN OCOD INTEGER NUMBER 
LESS THAN 32768 ON THE FIRST ENTRY TO GAUSS. THEFEAFT! 
IT WILL CONTAIN A UNIFORMLY OISTRIBLTED INTEGE? Faroa 
NUMBER GENERATED 3Y THE SUBROUTINE FOR USE ON THE NEY 
ENTRY TO THE SUBRUUTINE. 

S eTHE DESIRED STANDARD DEVIATION OF THE NORMAL 
OISTRIBUTICN. 

AM THE DESIRED MEAN OF THE NORMAL OISTRIBUTION 

Vo oeTHE VALUE OF THE COMPUTED NORMAL RANDOM VARTABLE 


REMARKS 
THIS SUBROUTINE USES RANOU WHICH IS MACHINE SPECIFIC 


SUBROUTINES ANO FUNCTION SUBPROGRAMS REQUIRED 
RANOU 


METHOO 
USES 12 UNIFORM PANDOM NUMBERS TO COMPUTE NORMAL RANDOM 
NUMBERS 8Y CENTRAL LIMIT THEOREM. THE RESULT IS THEN 
ADJUSTED TO MATCH THE GIVEN MEAN «0 STANDARD CEVIATION, 
THE UNTFORM RANOOM NUMBERS COMPUTED WITHIN THE SUBROUTINE 
ARE FOUND @Y THE POWER RESIDUE METHOD, 


SSO STHHHESFTSSHREEHREHHHESHSESSHEHEESEHSCEHEHHSHHC HCAS FEET EHFHoEHHSEHHEKB SEH EH HOS 


SUBROUTINE GAUSSCIX,S,A%,V) 
Aza, 

00 S@ Iste12 

CALL TANDUCIX,IY,Y) 

IxXslY 

AsAGY 

VE(A*6.9) *SF44 

RETURN 

END 


c-10 


NAAM 


4@ 


18@ 


1a@5 
110 


149 


SUBROUTINE GHK 


THIS SUBROUTINE COMPUTES THE VALUE OF Gy H, ANID K FROM BETA 


SUBROUTINE GHK (Gs He Ke Y) 
REAL Ke KY . 
INTEGER Ci 


THIS S3LOCK OETERMINES IF BETA wILL JE ENTERED FRO! THE KEYEQARR 
OR COMPUTED VIA A CALL TO SUAPOUTINE BETA, 


PRINT #, ° ENTER 1 IF YOU WANT TQ ENTER THE VALUE OF SETA FPR 
1 THE KEY8DARD* 
NB = @ 
C1 2 @ 
READ (5,44,ERR = 23) NB 
FORMAT (11) 
IF ( N@ .NE. 1) THEN 
CALL SETA (3) 
ELSE 
PRINT *,° PLEASE ENTER THE VALUE CF SETA’ 
READ (5,86,ERR = 62) 8B 
FORMAT (E16,9) 
ENDIF 


COMPUTE THE VALUE OF G, H, AND K, RETURN VALUE IN Y 


8 
1. 2 Y wx 
3.0/2. * ( € 1. © Yo) w (C{,°Y) ) 
1/25. * € CLee¥) #3) 
PRINT 190, Y 
: FORMAT (15X,° THE VALUE OF BETA IS °%,&15.9/) 
. WRITE (61298) Y 
PRINT 119, Ge He K 
FORMAT (° GIS °-E12.5,/° 4 IS °,E1L2.5,/ 
1° K IS %,€12.5//) 
PRINT *, ° ENTE? 1 TO CHANGE VALUF OF BETA’ 
acce’T 144, C1 
IF (C1 .EA, 1) GO TO 22 
FORMAT (11) 
RETURN 
END 


xXIO~< 
unin 


C-:1/(C-12 blank) 


DISTRIBUTION 
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